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Abstract 



Entropic entanglement measures of a two-dimensional system of two Coulombically interacting particles confined in an 
t-H anisotropic harmonic potential arc discussed in dependence on the anisotropy and the interaction strength. The harmonic 

approximation appears exact in the strong interaction limit, allowing determination of the asymptotic expression for the 

^ ^ 'linear entropy. Entanglement properties are dramatically influenced by the anisotropy of the confining potential in the 

&Q,strong-correlation regime. 
£3 . 
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1. Introduction 

The Hookean system composed of Coulombically in- 
teracting particles confined in a harmonic potential is of 
increasing interest as it provides an effective model of 
Semiconductor quantum dots (QDs) Determination 
of the amount of entanglement in various states of such 
systems is important in view of their possible applica- 
tions in quantum information technology Q. The sim- 
plest candidate for studying the entanglement properties 
is the two-particle Hookean system. Although the system 
was considered in various theoretical contexts both in 2D 
0, H S, @, 0, @, and in 3D case [1, E3, EI El , the influ- 
ence of the confinement anisotropy on entanglement has 
not been investigated so far. In this paper we undertake 
the investigation of this issue, restricting ourselves to the 
2D case, where the Hamiltonian is of the form 



Pi 

2m* 



I 2 1 



e*|r 2 - ri 



With e* being the effective dielectric constant and m* 
the effective electron mass the above Hamiltonian is a fre- 
quently used model of the two-electron QD. 

2fi_ E ^ hu^E the 
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After transformation r n- 
Schrodinger equation takes a form 

fftf(ri,r a ) =E*(n,r 2 ), 

where the Hamiltonian is given by 



H 



2 

£[ 
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2x1 + 2e 2 Z/?] 



r 2 - ri 



(3) 



The dimensionless coupling g = |y y -^j^ represents the 
ratio of the Coulomb repulsion to the confinement energy 



and the dimensionless parameter e = — measures the 
anisotropy of the confining potential. 

We will analyze the dependence of the entanglement be- 
tween the particles in the ground-state of the system on 
the interaction strength g and the anisotropy parameter e, 
paying particular attention to the regime of large g. In the 
case of finite e, the limit of g — > oo corresponds to the situ- 
ation in which both frequencies of the trap tend to zero. In 
this case, regardless of the value of e, the correlations play 
an essential role. In this paper we provide a method for 
determining the natural orbitals in the g — > oo limit by ap- 
plying the harmonic approximation to the anisotropic con- 
finement case. We derive an explicit representation of the 
asymptotic natural orbitals in terms of one-dimensional 
orbitals defined by integral equations. This enables easy 
determination of the asymptotic occupancies and entangle- 
ment entropies for the whole range of e. For finite values 
of g we determine the numerically exact results with the 
^Rayleigh-Ritz method and demonstrate how the asymp- 
totic values are attained. For all anisotropics, including 
the isotropic limit e = 1, the asymptotic values of entangle- 
ment entropies are properly determined by the harmonic 
approximation. 

The paper is arranged as follows. In scction[5]wc discuss 
the two-particle state characteristics. In section|3]wc show 
the reliability of the harmonic approximation in the regime 
of large g and provide the asymptotic Slater-Schmidt de- 
composition. In this section the entanglement properties 
are examined in detail. Finally, in section 0] we make the 
concluding remarks. 

2. Two-particle state characteristics 

2.1. Energy eigenspectrum 

Consider a 2D system consisting of two identical 
fermions with a Hamiltonian given by ([3]). Since the 
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Hamiltonian does not depend on spin, the solution of the 
Schrodinger equation 



ff¥(Ci,C2) = C2), 
where Q — (r^crj), factorizes to the form 

* ? (Ci,C 2 )-X?^ ± (ri,r 2 ), 



(4) 



(5) 



where s z = a\ + 02 and the labels correspond to the 
singlet (S) and triplet (T) states, the spin functions of 
which are given by xZ=o = 75 {\h>\~h> T I — i>l i>) and 

Xj =±1 = l ±, 3>l ± i>' The spatial wavefunctions ip^, that 
are symmetric (+) or antisymmetric (— ) under permuta- 
tion of the electrons, may be chosen real, since the inter- 
action and confinement potentials are real functions. 

Introducing the center of mass R = i(ri + r 2 ) = (X, Y) 
and relative coordinates r = r 2 ri = (x,y), the Hamil- 
tonian ((3]) may be written as H = H^ 1 + H r , where 



H 



H 



R 



-V 2 R /4 + 4(X 2 + e 2 Y 2 ^ 
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(6) 
(7) 



With wavefunction represented as a product 
V{x,y,X,Y) = ^ R (X,Y)tp r (x,y), the Schrodinger 
equation @ separates into two equations 



ffV(R) = fiV(R)> 

H T i/) r (r) = ^ r ^ r (r), 



(8) 
(9) 



where the total energy E = E r + E R . As the center-of- 
mass coordinate remains the same upon the interchange of 
electrons, the symmetry requirement reduces to the sym- 
metry of the relative wavefunction under inversion r > — r. 
Because of the invariance of H r to reflections about the 
x~ and y— axes, the (x, y)— parity of tp r (x,y) is well de- 
fined. The parity (+, +) or (— , — ) corresponds to spin 
singlet cigenfunctions, and the parity (+, — ) or (— ,+) to 
the spin triplet ones. The solutions of the CM equation 
(|8) are given by 



(10) 



i>* m {X,Y) = y n {2X)y m {2jiY), 

with the HO eigenfunction 

<p n (x) = e- x2 / 2 H n (x), (11) 

and corresponding energies E R n = 2(n + |) + 2e(m + |). 
The relative-motion equation ((5J) is separable only if e = 
1, 2, 1/2. The first case, oj x = u> y , being separable in polar 
coordinates, was much studied in the literature (13. HBj 
and in addition the closed-form solutions for particular 
frequencies have been derived flij ]. In the case ui y — 2ui x , 
Eq.© is separable in parabolic coordinates [13] and the 
closed-form solutions may be also obtained for particular 
values of to x Q • For other ratios of confinement frequencies 
the 2D Schrodinger equation has to be solved by numerical 
techniques. 



2.2. Entanglement measures 

Entanglement is a term used to describe quantum corre- 
lations between the particles. A convenient tool to analyze 
those correlations is the reduced density matrix (RDM) 
defined as [181 



p red (C,C)=Tr C2 ((C,C 2 |*)(*|C2,C})- 



(12) 



For identical fermions the bi- partite pure state Y$>) can 
be expressed as a combination of the Slater determinants 
made out of one-particle spin-orbitals in which its RDM 
IP|) is diagonal [H, 0, 



21]. The number of expansion 



coefficients appearing in the Slater decomposition, that 
are different from zero, is called the Slater rank (SR). The 
pure state of identical particles is considered nonentanglcd 
if SR = 1, i.e. the only correlations that exist between 
the fermions can be attributed to their indistinguishable 
nature 



21] 



The total wavefunction factorizes into spatial and spin 
components (J^Jl and the same holds for the reduced density 
matrix (HI]) 



pted,sSC,C') = Pred,sM °' ^Pred ( r . V ) » 



(13) 



where p re(i (r,r ) = J ^ ± (r, r 2 )V' ± (r2, r )dr 2 and the spin 

P artS Pred,B.=0 = Pted,s z =0 = ^(1/2,1/2), pj ed . Sz = 1 = 

diag(l,0) and P^erf,^=-i = diag(0, 1). 

The eigenvalue problem for the spatial part p red (i",r ) 
can be written in the form 



pf ed {v,v)vf{v)dv = A±«j t (r), 



(14) 



which determines the natural orbitals v l and the occu- 
pancies Xf. The families {w ; + (r)} and {vf(r)} form or- 
thonormal basis sets in the space of symmetric and an- 
tisymmetric functions, respectively. The spatial parts of 
two-particle functions may be represented in terms of nat- 
ural orbitals. Representation for the symmetric function 
takes the form of the Schmidt decomposition 



■0 + (ri,r 2 ) = ^fc+i>+(ri)v+(r 2 ), 



(15) 



and in the case of antisymmetric function the Slater de- 
composition holds [2l[ 

V>~(ri,r 2 ) = ^^[%-i( r i) u 2i( r 2)-%-i( r 2)%( r l)]-( 16 ) 

The coefficients and k~[ are related to the eigenval- 



ues of (HU) by A+ = [k+] 2 and A" - 



respectively. 



The eigenvectors v 2l and v 2l _ 1 correspond to the same 
eigenvalue A ; ~ , which means that eigenvalues of the spatial 
RDM of antisymmetric wavefunction are doubly degener- 
ate and satisfy the conservation of probability 2 A^ = 
1. We assume that the eigenvalues A;* are ordered such 
that Aq > \f > .... It can be easily inferred from the ap- 
propriate decompositions of the spatial parts (fT5j) or (fTI 



2 



that the SR is related to the number n of non-vanishing 
eigenvalues of the spatial RDM (being not necessarily dif- 
ferent) as follows: the singlet and the triplet states with 
s z = have SR = n, whereas the triplet states with 
s z = — 1,1 have SR = n/2. The two-particle state is non- 
entangled if and only if SR = 1 and deviations from such a 
form may be used to measure the amount of entanglement 
in the system. 

The entanglement depends on the whole spectrum of the 
RDM and its characteristics may be constructed from the 
natural orbitals occupancies A; 



241. In this work 



we consider the von Neumann (vN) entropy and the linear 
entropy that are the most popular measures of entangle- 
ment in pure states. The vN entropy [25, 26[ 21 1 is defined 

as 

S = -Tr[p red Log 2 Pred}, (17) 

and in the two-particle case separates into 

S = S spm + S vN , (18) 

where the spin contribution S spm = 1 if s z = and 
gspin = o if Sz = ±L 

The space part that depends on interactions may be 
calculated with the RDM eigenvalues as 



^ XiLog 2 \i 



1=0 

The linear entropy [2] 



L = 1 - Trp 



red 



(19) 



(20) 



can be calculated directly from the above definition, which 
is a big advantage over the case of vN entropy (JTTJ) since 
diagonalization of the RDM is not needed. 



3. Harmonic approximation 

3.1. Approximate spectrum 

The harmonic approximation was applied with success 
in the case of circular [H, 28, 2|| and anisotropic confine- 
ment potential [30j ]. The anisotropic potential V(x,y) = 



y 2 has two local minima at (x c i , 0) 
and (— x c i,0), where x c i = (g/2)i is the classical equilib- 
rium distance between the Coulombically interacting par- 
ticles in the trap. Expanding V(x, y) into a Taylor series 
around the minimum at r m ,„ = (x c i , 0) and retaining the 
terms up to second order, the relative motion equation @ 
gets approximated by the Schrodinger equation 



A r +V(r mm ) + 3{x - x cl f + (e 2 



l)y 2 W 



the eigenvalues of which are of the form 

K m = V{v mm ) + 2^{\+n) + 2{\ + m)^e 2 ~l.{22) 



The corresponding cigenfunctions 

r nm (x, y) = <p n (3 1/4 (* ~ Zc/))Vm((e 2 - l) 1/4 y), (23) 

where tp n is given by (|11[) . provide approximations to 
the relative motion wavefunctions only around (x c ;,0). 
The approximations around (~x c i,0) are obtained by the 
transformation x i— > —x. The states of relative motion 
with well-defined parity (x, y) may be constructed as 



(24) 



where the sign +/— relates to the even/odd x— parity. No- 
tice that each energy level of the relative motion in the 
harmonic approximation is doubly degenerate, namely the 
states with spatial x-y parity: (— , +), (+, — ) arc degener- 
ate with those of (+, +), (—,—), respectively. 

The quality of harmonic approximation has been tested 
by comparison with numerical solutions of the relative mo- 
tion equation ^ obtained through the exact diagonaliza- 
tion method in the basis of the two-dimensional harmonic 
oscillator eigenfunctions of the required parity. In Fig. [T] 
the excitation energies E\ — Eq of the Hamiltonian (O are 
plotted in function of In g for a large (e = 1.7) and a small 
(e = 1.01) anisotropy. The horizontal lines mark the re- 
sults of the harmonic approximation (|22[) . which becomes 
exact at g — y oo. We observe that for smaller values of e 
the asymptotic behavior of the relative motion energies is 
reached at larger values of g. This is consistent with the 
fact that the closer the value of e to one, the larger the 
value of g at which the harmonic approximation is appli- 
cable. 

For transparency of presentation we concentrate on an- 
alyzing the singlet and triplet states of the lowest energy 
that become degenerate in the limit of g — > oo. They cor- 
respond to the ground state of CM; thus the total spatial 
wavefunctions are ip + = V'So^oo an< ^ V' - = - 0ob~' ( /'oo' 



re- 



spectively Harmonic approximations to those functions 
can be written in the form convenient for further analysis 



ip ± (xi,x 2 ,yi,y2) = C±(g,e)h(yi,y 2 )(q(xi,x 2 ) ± q(x 2 ,xi)), (21 



where 



q(xi,x 2 ) = ip (3 1/4 (x 2 -x±- x c i))ipo(x 2 + xi) = 

e -^lV3(x 2 -xi-x cl ) 2 + (ani+a:2) 2 ] ^ ^6) 



h(yi,y 2 ) = <A)((e 2 ~ l) 1/4 (2/2 - 2/i))<A)(Ve(2/2 +Vi)) = 



,(27) 



(21) and the normalization constants are given by 



V2W(6 2 -1)^ 
C±{g,e) = ; g 

1±e -V3(f)s 



,v/23M(e 2 -l)s 



(28) 
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Figure 1 : Excitation energies Ef — Eq in function of In g for e = 
1.7 (top) and e = 1.01 (bottom). The horizontal lines mark the 
asymptotic results of the harmonic approximation. 

3.2. The linear entropy 

Quantum entanglement in two-particle systems is usu- 
ally quantified by entropic measures. The easiest to cal- 
culate is the linear entropy that may be obtained directly 
from the definition (|20p using the RDM calculated from 
the spatial wavefunction. The asymptotic expression for 
wavefunction (|25[) enables us to determine the linear en- 
tropy in the limit of g — > oo. For the lowest singlet state 
the explicit expression reads 

3*(e 2 -l)We(l-^) 

L»-~(e) = 1 / (29) 

e + v e — 1 

Although the above formula is strictly valid only if e > 1, 
since it has been derived from the harmonic approxima- 
tion, we observe that its limit at e — > 1 + is equal to one 
and coincides thus with the exact asymptotic value of the 
linear entropy. 

In FigJS] we show how the linear entropy approaches 
the asymptotic limits for various values of the anisotropy 
parameter. The results at finite values of g have been 
obtained by numerical integration with the use of wave- 
functions determined by the Rayleigh-Ritz procedure. At 



large anisotropy the dependence on g is very similar and 
the asymptotic values of the linear entropy for e = 2 and 
e = 4 are nearly equal to that in the infinite anisotropy 
limit L 9_Kx> (e -> oo) « 0.759142. Only at e < 1.4 asymp- 
totic values of the linear entropy start to visibly differ from 
0.759142 and reach the value one at e — > 1 + . 

In this range of e the linear entropy exhibits a local 
maximum and approaches monotonically from above its 
asymptotic value determined by the formula (|29[) . which 
confirms its validity. The smaller is the anisotropy, the 
larger is the value of g at which it occurs, which is in 
accordance with our earlier discussion. The dependence of 
the linear entropy on e becomes strong for large enough g, 
but only in a small vicinity of e = 1. 




Lng 

Figure 2: The linear entropy of the lowest singlet state for various 
anisotropy ratios e in function of lng. The horizontal lines show the 
asymptotic values at g — ¥ oo. 



3.3. Asymptotic Slater- Schmidt decomposition 

In order to calculate the vN entropy, determination of 
the natural orbital occupancies is necessary. The asymp- 
totic behavior at g — » oo may be determined by deriving 
the Schmidt decomposition of the asymptotic two-particle 
function (|25p. To this end it is convenient to introduce 
new coordinates x\ x\ — x 2 t— > i,i + ^i, so as the 
function q (|26[) transforms into 

q(x 1 ,X 2 ) ^ q(x 1} X3) = e -5[v / 3(i2-ii) 2 +(£i+i2) 2 ]^ ( 3Q ) 

Being real and symmetric, the function q has the Schmidt 
decomposition 

q( Xl ,x 2 ) = ]T k^^\ Xl )^\x 2 ). (31) 

n=0 

It is worthwhile to notice that the function (f3"U)) is g- 
independent and so are the orbitals &n Emd the corre- 
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sponding expansion coefficients kn\ They can be deter- 

(32) 



mined from the integral equation 

q(x,X )^ ) {x)dx =kP#W(x), 



Changing the variables back to x\ and X2, one gets 

q( Xl ,x 2 ) = ^^^'(^i + ¥)€H*2 " (33) 



The function /i(yi,y 2 ) is real and symmetric, thus its 
Schmidt decomposition reads 



m=0 



(2) (2) 

where i) m and fcm are determined by 



h{y,y')0g\y')dv =k$0W(y), 



(34) 



(35) 



(2) 

with the orthogonal orbitals fim assumed to be normalized 
to unity. Using the expansions (|3"5f and (fMf we represent 
the wavefunctions ([23|l as 

^ ± (n , r 2 ) = £ fc± [L, (n (r 2 ) ± R t (r 1 )L t (r 2 )] , (36) 



2=0 



where Z = (n, m), fc* = C±(g, e)fc„ fcm and 

L i (r) = ,? i ^ + -|i )y )=^)(x+^)^)( 2/ ) ) 

^(r) = ^-^ L ,y) = ^ 1) («-^)^ ) (y). 

are the one-particle orbitals centered around classical equi- 
librium points which satisfy (Li\Lk) = Ski and = 
Ski ■ If we define the new orbitals 

v l (r)= Rl{r)+ r - Ll{T) Mr)= L ^- R ^ 



V2 



V2 



that fulfil (ui\vi) = 0, the spatial wavefunctions can be 
expressed as 

^ + (ri,r 2 ) = ^2 k+[vi(Y 1 )vi(r 2 ) - uj(ri)itj(r 2 )], (37) 



2=0 



and 



V> (n,r 2 ) = ^);, [tii(ri)B|(r2)-»i(ri)u ( (r 2 )]. (38) 

2=0 

In the limit of g — > oo we have fc ; + = = ki and 
the integral overlap (Li\Rk) vanishes for any I, k, and be- 
cause of that HV^H 2 = ^Y.i k h ( u i\ u k) = (vi\v k ) = 5ik, 
(ui\vk) = 0. In this limit the symmetric spatial wave- 
function coincides with the absolute value of the antisym- 
metric one ■0 + (ri,r 2 ) = |^~(ri,r 2 )|. Moreover, (|37p and 
(|38l) yield the Schmidt decomposition (|T5"T) and the Slater 
decomposition (|16D , respectively. In both cases the decom- 
position of the spatial RDM reads 



P 



£ d (r,r') =J2kf[vi(r)v l (r')+ui(r)ui(v)], (39) 



which shows that the eigenvectors vi and ui correspond to 



the same occupancy Aj 



23g£?(e 2_ 1)g (1) (2) 



{kX'k^y with the 



coefficients kn and kn 1 determined by Eqs. (|32p and ([55)1 . 
respectively. The above equations can easily be solved 
through discretization technique. By discretizing the x(y) 
and x (y ) variables with equal subintervals of length A.t 
(Ay), the integral equations turn into algebraic eigenvalue 
problems 



,(2) 



E 



[A^-S rs k^(z s ) = 0, r,i = 1,2, 



(40) 



where A)-s = 9(^V,^ s )Aa; and A^s =- h(y r ,y s )Ay. Di- 
agonalization of the matrix [^4^2 ] n x at provides thus a set 
of approximations to TV of the largest modulus eigenval- 
ues k„ . One has to notice that the function h depends 
on e and the range of argument where its value is non- 
negligible strongly increases when e — > 1 + . Therefore, in 
the case of small anisotropy, an appropriately large inter- 
val must be discretized to achieve a reasonable accuracy 

( 2) 

in the diagonalization of [Ar S ]NxN- In Fig- [3] the four 
largest natural orbital occupancies determined in function 
of g by the numerically exact Rayleigh-Ritz procedure are 
shown for the nearly isotropic example of e = 1.01. It is 
seen how they converge into asymptotic doublets of the 
values determined in the harmonic approximation by the 
procedure described above. The behavior of the asymp- 




Figure 3: Four largest occupancies in function of In g at fixed 
anisotropy e = 1.01. The asymptotic values are marked by hori- 
zontal lines. 

totic occupancies at varying anisotropy ratio e is shown 
in Fig. 21 For a better display the results are presented 
in function of ln(e — 1), since the occupancies are strongly 
sensitive to small changes in e around e = 1. In the limit 
of ln(e — 1) — > — oo, the values of all asymptotic occupan- 
cies tend to zero in such a way that their sum is equal 
to i. Notice that this behavior is related to the degener- 
acy that appears in the energy spectrum of the system in 
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the circular symmetry limit. The clustering of occupancies 
around zero value is clearly visible in Fig|3J With increas- 
ing anisotropy the situation changes dramatically, namely 
a tiny deviation from circular symmetry leads to the lift- 
ing of the clustering. The value of ln(e — 1) s» —10, at 
which the effect is clearly visible, corresponds to a very 
small anisotropy, e w 1.00005. At higher confinement 
anisotropics, all the occupancies but the largest one exhibit 
a local maximum. The larger the value of the occupancy, 
the larger is the value of e at which the maximum occurs. 
Above the critical threshold e cr « 1.4 (ln(e cr — 1) « —1) 
all those occupancies saturate at vanishingly small val- 
ues. The largest occupancy Aq - * 00 performs differently, 
as it grows monotonically with increasing e at the cost of 
the remaining occupancies and saturates above e cr . The 
limit of ln(e — 1) — > oo corresponds to e — > oo, i.e. a 
one dimensional motion along the x-axis. In this limit 
Aq - * 00 « 0.490688 and the sum of all the remaining occu- 
pancies 2 J2i=i Af^ 00 is only about 0.018624, which means 
that the natural orbitals vq(t) , uq(t) are the only two that 
are substantially occupied. This indicates that in the case 
of strong interaction for enough anisotropic confinement 
(e > e cr ), the terms higher than I — in (j36|) contribute 
very little and the spatial functions approach the form 

ift^tri.ra) « -^[L (ri) J Ro(r 2 )±i?o(ri)io(r 2 )].(41) 

The total singlet and triplet wavefunctions with s z = are 
thus well approximated by 

*£=o(Ci,G0 ~ ^(\i>\-i>T\-i>\i>) x 

(L (ri)i2o(r 2 ) ± i2o(ri)L (r 2 )). (42) 

Each of them constitutes a sum of two Slater determinants 
(SR = 2) and represents an entangled state. The situation 
is different for the triplet components with s z = ±1, since 
their total wavefunctions 

C= ± i(Ci,C 2 )« 

-^=|±|>|±|>(io(ri)i?o(r 2 ) - i? (ri)i (r 2 )) (43) 

represent one Slater determinant (SR = 1) and those 
states have to be regarded as non-entangled. It has to 
be stressed however, that this is only an approximate re- 
sult since even at e — s- oo the other occupancies do not 
strictly vanish. 

3.4- Comparison of the linear and vN entropy 

Since the linear entropy is relatively easy to calculate, 
several works proposed jij, [26[ |3l|, [32| to use it instead of 
the vN entropy arguing that their dependence on parame- 
ters of the system is similar. Here we test this conjecture 
for the lowest singlet state by comparing the behavior of 
both measures in function of e at g — > oo, where the dif- 
ferences between both measures are the most pronounced. 
We observed that for e large, the dependence between the 
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Figure 4: Behavior of the six largest asymptotic occupancies in func- 
tion of ln(e— 1). The limit of the largest occupancy Ag^°° « 0.490688 
as e — ► oo is marked by a horizontal line. 

vN and linear entropy is approximately linear and well fit- 
ted by S 9 ^°°(e) w 12.8L 9 ^°°(e) - 7.6. The comparison of 
the so rescaled linear entropy with vN entropy is presented 
in Fig. [5] in function of ln(e— 1). Above e cr w 1.4 both en- 
tropies saturate at constant values, Sf~^ « 2.13618 and 
^e^too ~ 0.75914, respectively. The closeness of those val- 
ues to 2 and 3/4 can be attributed to the dominance of the 
two occupancies A^ ps A^ in this regime. The vN entropy 
grows with decreasing e, which reflects the fact that the 
smaller is the value of e, the larger number of orbitals in 
the sum (|36p becomes important. The increase is nearly 
linear which means that the asymptotic vN entropy varies 
logarithmically with e and tends to infinity when the con- 
finement becomes circularly symmetric (e — > 1 + ). Being 
a bounded function, the linear entropy performs differ- 
ently at very small anisotropics, where its behavior starts 
to deviate from the linear one and approaches the maxi- 
mum possible value 1 in the limit of e — > 1 + . Although 
in different ways, but both entanglement measures clearly 
demonstrate that the more circular is the confinement of 
strongly interacting particles, the more entangled is the 
system. The difference in behavior of the vN entropy and 
the linear entropy appears only in the vicinity of e = 1, 
where the system becomes degenerate. 

4. Summary 

We performed a detailed examination of the entangle- 
ment properties for the system of two Columbically inter- 
acting electrons confined in a 2D anisotropic harmonic po- 
tential. The harmonic approximation has been developed 
in order to study the strongly interacting (weak confine- 
ment) case. Using the harmonic approximation we derived 
an explicit expression for the linear entropy of the lowest 
singlet state at g — > oo. The occupancies A^ 00 in the 
asymptotically degenerate lowest singlet and triplet states 
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12.8L^°°-7. 




Figure 5: Comparison of the asymptotic vN entropy S 9— > °° and the 
(appropriately rescaled) linear entropy L 9 ~*°° in function of ln(e— 1). 



may be easily determined numerically within this approx- 
imation. Performing numerical calculation through the 
Rayleigh-Ritz method we have also calculated the charac- 
teristics of the system at finite values of g. We demon- 
strated how the occupancies reach their asymptotic values 
as g is increased. The asymptotic occupancies are strongly 
sensitive to changes of e in the range of 1 < e < e cr f=s 1.4, 
while above e cr they practically reach the values corre- 
sponding to infinite anisotropy e — > oo, where the value of 
Ag - * 00 w 0.490688 is the only substantial. This results in 
non-trivial entanglement properties since only two natural 
orbitals contribute significantly to the lowest singlet and 
triplet states. Furthermore, we have verified that the lin- 
ear entropy is almost linearly related to the corresponding 
vN entropy except of the region of very small anisotropy 
e < 1.00005. Our calculations have shown that the entan- 
glement is relatively insensitive to the shape of the har- 
monic confinement if the interaction is very weak (strong 
confinement case). In the strong correlation case the in- 
fluence of anisotropy is much more pronounced. 

It would be interesting to carry out an analysis of en- 
tanglement for higher excited states. The thermal entan- 
glement and the interaction with the environment will be 
the topic of further studies. 
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